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Pseudoscalar Meson Radial Excitations 
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Goldstone modes are the only pseudoscalar mesons to possess a nonzero leptonic decay constant 
in the chiral limit when chiral symmetry is dynamically broken. The decay constants of their 
radial excitations vanish. These features and aspects of their impact on the meson spectrum are 
illustrated using a manifestly covariant and symmetry preserving model of the kernels in the gap 
and Bethe-Salpeter equations. 
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The meson spectrum contains three pseudoscalars 
[I G (J P )L = l-(0-)S] with masses below 2GeV Q: 
tt(140); tt(1300); and tt(1800). Of these, the pion [tt(140)] 
is well known and much studied. The other two are ob- 
served, e.g., as resonances in the coherent production of 
three pion final states via pion- nucleus collisions 2J. In 
the context of a model constituent-quark Hamiltonian, 
these mesons are often viewed as the first three members 
of a QQ n ^o trajectory, where n is the principal quan- 
tum number; i.e., the 7r(140) is the S-w&ve ground state 
and the others are its first two radial excitations. By 
this reasoning the properties of the 7r(1300) and 7r(1800) 
are likely to be sensitive to details of the long range part 
of the quark-quark interaction because the constituent- 
quark wave functions will possess material support at 
large interquark separation. Hence the development of an 
understanding of their properties may provide informa- 
tion about light-quark confinement, which complements 
that obtained via angular momentum excitations |3j. 

This view might reasonably be held about all trajec- 
tories of radial excitations. However, the pseudoscalar 
trajectory is particularly interesting because its lowest 
mass member is QCD's Goldstone mode. An explana- 
tion should therefore simultaneously describe: (1) chiral 
symmetry and its dynamical breaking; and (2) the pos- 
sible correlation of the trajectory's higher mass members 
via an approximately linear radial Regge trajectory |_l]. |o| . 
Outcome (2) does not require that confinement in light- 
quark systems be expressed through the formation of a 
flux tube |6(. It is easily obtained in Poincare invariant 
quantum mechanics pj but requirement (1) is not. 

A Poincare covariant and symmetry preserving treat- 
ment of quark-antiquark bound states can be based on 
the Bethe-Salpeter equation (BSE) 



r t u(k;P)= / [ X (q;P)]srK u s (q,k;P), 



(1) 



where: k is the relative and P the total momentum of the 
constituents; r,. . . , u represent colour, Dirac and flavour 
indices; xfe P) ■= S(q + )r(q;P)S(q.), q± = q ± P/2; 
and J represents a translationally invariant regulari- 
sation of the integral, with A the regularisation mass- 
scale |a ■ I n Eq. , S is the renormalised dressed- 



quark propagator and K is the fully amputated dressed- 
quark-antiquark scattering kernel. The product SSK 
is a renormalisation point invariant. Hence, when the 
kernel is expressed completely in terms of renormalised 
Schwinger functions, the BSE's solution is independent 
of the regularisation mass-scale, which may be removed; 
viz., A — > oo. 

In a given channel the homogeneous BSE only has solu- 
tions for particular, separated values of P 2 : P 2 = —m^, 
where m n is a bound state's mass, whereat T n (k;P) is 
that bound state's Bethe-Salpeter amplitude. (We use a 
Euclidean metric, with: {7^,7,,} = 25^; 7^ = 7^; and 

a- b = J2i=i a ibi- For a timelike vector P^, P 2 < 0.) 
In the flavour nonsinglet pseudoscalar channel the lowest 
mass solution is associated with the 7r(140). The homo- 
geneous BSE next possesses a solution when P 2 assumes 
the value associated with the mass of the 7r(1300). This 
pattern continues so that in principle one may obtain the 
mass and amplitude of every pseudoscalar meson from 
Eq. 0. Herein we will illustrate this in practice for the 
two lowest-mass flavour-nonsinglet pseudoscalar mesons. 
The dressed-quark propagator appearing in the BSE's 
kernel is determined by the renormalised gap equation 



E(p), 



(2) 



Sip)- 1 = Z 2 (i~f P + m b > 

S(p) = Z 1 J 9 2 D^(p~q)^ ll2 S(q)T:( q ,p),(3) 

where D^ v is the dressed-gluon propagator, T u (q,p) 
is the dressed-quark-gluon vertex, and m bm is the A- 
dependent current-quark bare mass. The quark-gluon- 
vertex and quark wave function renormalisation con- 
stants, £/i,2(C 2 i A 2 ), depend on the renormalisation point, 
£, the regularisation mass-scale and the gauge parameter. 
The gap equation's solution has the form 

(4) 

It is obtained from Eq. J2J augmented by the renormali- 
sation condition 



Sip)- 1 = i 1 -pA{p 2 ,( 2 )+B{p 2 X 2 ). 



Sivy 1 



P 2 =C 2 



ij-p + miC), 



where m(£) is the renormalised mass: 

Z 2 (C 2 ,A 2 )m bm (A) = Z 4 (C 2 ,A 2 )m(C), 



(5) 



(6) 



with Z^ the Lagrangian mass renormalisation constant. 
In QCD the chiral limit is unambiguously defined by 

Z 2 (C 2 ,A 2 )m bm (A)=0,VA»C, (7) 

which is equivalent to stating that the renormalisation- 
point-invariant current-quark mass rh = 0. 

The 1 _ (0 _ )5 (n 1 ^) trajectory contains the pion, 
whose properties are fundamentally governed by the 
phenomenon of dynamical chiral symmetry breaking 
(DCSB). One expression of the chiral properties of QCD 
is the axial-vector Ward-Takahashi identity 



P M rUfc;P) = S~ 1 {k + )i lh - 
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-2tm(0n(k;P), (8) 

which we have here written for two quark flavours, each 
with the same current-quark mass: {r l : i = 1,2, 3} are 
flavour Pauli matrices. In Eq. IJHJ), T^ (k; P) is the axial- 
vector vertex: 






(k;P) 



757m" 



+ / lxUq;P)lsrK™(q,k;P),(9) 



and r J 5 (k; P) is the pseudoscalar vertex 



[r 5 (*;;P)L 



75y 



{xUq;P)}srK™( q ,k;P). (10) 



The quark propagator, axial-vector and pseudoscalar 
vertices are all expressed via integral equations; i.e., 
Dyson-Schwinger equations (DSEs). Equation JSJ) is an 
exact statement about chiral symmetry and the pattern 
in which it is broken. Hence it must always be satisfied. 
Since that cannot be achieved veraciously through fine 
tuning, the distinct kernels of Eqs. (0), ©, ©, {TOl must 
be intimately related. Any theoretical tool employed in 
calculating properties of the pseudoscalar and pseudovec- 
tor channels must preserve that relationship if the results 
are to be both quantitatively and qualitatively reliable. 

A weak coupling expansion of the DSEs yields pertur- 
bation theory and satisfies this constraint. However, that 
truncation scheme is not useful in the study of bound 
states nor other intrinsically nonperturbative phenom- 
ena. Fortunately at least one nonperturbative systematic 
and symmetry preserving scheme exists. (Reference [lOJ 
gives details). This entails that the full implications of 
Eq. JSJl can be both elucidated and illustrated. 

Every flavour nonsinglet pseudoscalar meson is exhib- 
ited as a pole contribution to the axial-vector and pseu- 
doscalar vertices M: viz., 



n,(k;P) 
%Ti(k;P) 



- p2 + m J„~o 

P 2 +ml «0 



^t_TiJk-P), (11) 



P 2 + m 

Ptt 

P 2 + 



-rUkiP), (12) 



wherein we have omitted terms regular in the neighbour- 
hood of the pole, r£ (k; P) is the bound state's canoni- 
cally normalised Bethe-Salpeter amplitude: 

Ti n (k; P) = ^75 [iE nn (k; P) + 7 ■ PF*„(k; P) 
+ i-kk-PG« n (k',P) + <r lu ,k lt P v H Kn (k',P)]\(13) 

and 



U n 6 ij P» = Z 2 trj ymvXijQlP), (14) 
ip, n (0S iJ = Z 4 tr ^ l5X i(q;P). (15) 



The residues expressed exactly in Eqs. ifHH) and (fL5|) . are 
gauge invariant and cutoff independent. 

For a structureless pseudoscalar meson, F 7!n (k; P) = 
= G ff „(fc;P) = H Wn (k;P) in Eq. (JSJ. The" first two 
of these functions can be described as characterising 
the pseudoscalar meson's pseudovector components; and 
the last, its pseudotensor component. The associated 
Dirac structures necessarily occur in a Poincare covari- 
ant bound state description and they signal the presence 
of quark orbital angular momentum. 

Equation © combined with (JTTJl - (JTHJl yields H, 



/*•„«•** 



2m(C)ftr B (C) 



(16) 



i.e., an identity valid for every - meson and irrespective 
of the magnitude of the current-quark mass |ll| . 

We focus now on the ground state n = pseudoscalar 
meson and assume that all pseudoscalar excitations are 
more massive. In this case DCSB in QCD entails Q , via 
Eqs. I©,© and JT3J - (TSJ, 

f° o E„ (k;0)=B(k 2 ,( 2 ), (17) 

where /° := lim^— >o / 7rn , from which follows 



Pl (C):-- 
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m—>0 
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■<«)?, 



wherein 



-{qq)l= lim Z 4 (C 2 , A 2 ) iV c tr D 



s°(?,C), 



(18) 



(19) 



is the vacuum quark condensate [12J. The scalar func- 
tions in Eq. J3J that characterise the dressed-quark prop- 
agator are both positive definite at spacelike momenta 
[131 lll| . This fact along with their ultraviolet behaviour 
in the chiral limit [l5j guarantee < —{q~q)9 < oo. Fur- 
thermore, Eq. IJ17II and related quark-level Goldberger- 
Treiman relations [8| involving A(k 2 ;^ 2 ) ensure that 
< f® < oo. Hence m^ — in the chiral limit, and in 
the neighbourhood of this limit the so-called Gell-Mann- 
Oakes-Renner relation is a corollary of Eq. i|l(j|) . 

We next consider the n > pseudoscalar mesons: 
TO 7r„>o > rn^g by assumption, and hence m nn>0 ^ in 



the chiral limit. The existence of a bound state entails 
that x„(fc;P) is a finite matrix-valued function. More- 
over, the ultraviolet behaviour of the quark- antiquark 
scattering kernel in QCD guarantees that Eq. (|15|l is cut- 
off independent. Thus 



pUO 



lim pn n (C) < oo , V n . 

m— >0 



(20) 



Hence, it is a necessary consequence of chiral symmetry 
and its dynamical breaking in QCD; viz., Eq. (|16(l . that 



/° =0,Vn>l. 



(21) 



This result is consistent with Refs. 16], as appreciated in 
Ref. H3. 

This argument is legitimate in any theory that has a 
valid chiral limit. It is logically possible that such a the- 
ory does not exhibit DCSB; i.e., realises chiral symmetry 
in the Wigner mode. This is the case, e.g., in QCD above 
the critical temperature for chiral symmetry restoration 
[l8| . Equation l|16[) is still valid in the Wigner phase. 
However, its implications are different. 

Without DCSB; namely, in the Wigner phase, one has 



B w (0,( 2 )<xm{()(xm; 



(22) 



i.e., the constituent-quark mass vanishes in the chiral 
limit. This result is accessible in perturbation theory. 
Equation (|17f) applies if there is a massless bound state 
in the chiral limit. Suppose such a bound state persists 
in the absence of DCSB. (If that is false then considering 
this particular case is unnecessary. However, it is true at 
the transition temperature in QCD [T3.) It then follows 
from Eqs. JE! and J23 that 



'7TO 



-TV. 



(23) 



In this case the leptonic decay constant of the ground 
state pseudoscalar meson also vanishes in the chiral limit. 
It is always true that 



UoP*o(0 « 



■(&)$ 



(24) 



In the Wigner phase (12|, (qq)® w oc to 3 . Hence, via 
Eq. I|16fl . if a rigorously chirally symmetric theory pos- 
sesses a massless pseudoscalar bound state then 



,w 



(25) 



In this case there is also a degenerate scalar meson part- 
ner whose mass behaves in the same manner. 

We have elucidated exact results. They can be il- 
lustrated accurately in a model that both preserves 
QCD's ultraviolet properties and exhibits DCSB. For this 
purpose we employ the renormalisation-group-improvcd 
(RGI) rainbow-ladder DSE model introduced in Ref. Hal 
whose widespread application is reviewed in Ref. [20j. 



The heart of the model is an Ansatz for the Bcthe- 
Salpeter kernel in Eq. Q: 



K t r u s (q,k;P) = -g((k-q) 2 ) 



X Df v c {k q) 



7 "T 



7 "T 



(26) 



wherein D f ™ c (l) is the free gauge boson propagator and 



m = jL Dse -^ 



+- 



27T7„ 



111 



r+(l+s/A 2 QCD y 



F{8) 



(27) 



with T(s) = [1 - exp(-s/[4mf])]/s, m t = 0.5 GeV, r = 
e 2 - 1, 7 m = 12/25 and A QCD = A^ = 0.234. 

This form expresses the interaction as a sum of two 
terms. The second ensures that perturbative behaviour 
is correctly realised at short range; namely, as written, 
for (k — q) 2 ~ k 2 ~ q 2 > 1 — 2 GeV 2 , K is precisely as 
prescribed by QCD. On the other hand, the first term 
in G(k 2 ) is a model for the long-range behaviour of the 
interaction. It is a finite width representation of the form 
introduced in Ref. |2l| , which has been rendered as an in- 
tegrable regularisation of 1/fc 4 [22|. This interpretation, 
when combined with the result that in a heavy-quark- 
heavy-antiquark BSE the RGI ladder truncation is exact 
[lCj . is consistent with Q(k 2 ) leading to a Richardson-like 
potential [23| between static sources. 

The active parameters in Eq. I|27|l are D and u, which 
together determine the integrated infrared strength of 
the rainbow-ladder kernel, but they are not indepen- 
dent. In fitting a selection of observables (l9|. a change 
in one is compensated by altering the other; e.g., on 
the domain ui € [0.3, 0.5] GeV, the fitted observables 
are approximately constant along the trajectory toD = 
(0.72 GeV) 3 =: to 3 . Herein we use w = 0.38 GeV. 

With a given truncation of the BSE's kernel there is 
a unique kernel for the gap equation which ensures the 
axial- vector Ward-Takahashi identity, Eq. (J8J, is auto- 
matically satisfied 10] . This partner to Eq. (|26|) is a 
rainbow gap equation; i.e., Eq. J2J with 

^(p) = fg((p-qf)D^(p^q)^ ltl S(q)^. (28) 

We have calculated properties of the ground and first 
excited states. The first step was to solve the gap equa- 
tion, Eq. (|28|l , for a specified renormalised current-quark 
mass. With the dressed-quark propagator thus obtained, 
solutions of the homogeneous BSE, Eq. (JIJ with 126fl . 
were obtained via a straightforward numerical procedure 
that yields every scalar function in Eq. i|T3|l . The general 
procedure is described in detail in Ref. [3 , and Ref. J24J 
explains the method necessary to isolate excited states. 

In Fig. n we depict the lowest Chebyshev moments of 
the pseudoscalar amplitude in Eq. I|13|) : i.e., °' 2 E^ (k 2 ) 
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FIG. 1: Dimensionless low-order Chebyshev moments of the 
scalar function that characterises the dominant amplitude in 
Eq. 113H . for the ground (no) and first excited (ni) states. 



FIG. 3: Evolution of the mesons' leptonic decay constants 
with renormalised current-quark mass, m q = m(£ = 1 GeV). 
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FIG. 2: Evolution of the mesons' masses with renormalised 
current-quark mass, m q = m(C = 1 GeV). 



and 0,2,i E ni (k 2 ), where 



E^ 01 (k 
2 



dp sin 2 /3 Ui(cos (1) E W01 (k 2 ,k ■ P;P 2 ) . (29) 



with Ui(x) a Chebyshev polynomial of the second kind 
and k ■ P := cos /3V 'k 2 P 2 . The odd moments, i = 
1,3,5,..., etc., vanish in the case of equal mass con- 
stituents. The Chebyshev moments are obtained from 
the canonically normalised Bethe-Salpeter amplitudes 
but, for illustrative simplicity, the functions depicted are 
rescaled by the positive constant °E no (k 2 = 0) . 

All Chebyshev moments of E 7ri possess a single zero, 
whereas those of E„ exhibit none. This similarity to the 



wave functions of radial excitations in quantum mechan- 
ics is not particular to manifestly covariant BSE studies 
|25l l2r1 |. It is evident that the zeroth Chebyshev mo- 
ment almost completely determines the pseudoscalar am- 
plitude in the ground state pseudoscalar meson. For the 
first excited state, however, the second moment is also re- 
quired to obtain a good approximation to E 1Tl (k; P). The 
pseudovector and pseudotensor amplitudes are nonzero 
in the ground and first excited states, and they are ma- 
terially important in the calculation of their properties. 
The bulk qualitative features of the scalar functions char- 
acterising these amplitudes are the same as those de- 
scribed in connection with the pseudoscalar amplitude. 

In Fig.[21 we depict the evolution with current-quark 
mass of the masses of the ground and first excited pseu- 
doscalar states, and in Fig.Elwc display the behaviour 
of the leptonic decay constants. Calculated results at 
points of particular interest on these trajectories are pre- 
sented in Table [I] NB. The ladder truncation supports 
ideal mixing. This is not a good picture of ground 
state ss pseudoscalars, for which s-channel gluon (OZI 
suppressed) contributions are important. (See, e.g., 
Ref. [23 •) However, it should prove increasingly reliable 
as the current-quark mass increases, since the RGI ladder 
truncation is exact in the static source limit, and/or as 
the mass of the bound state increases because this mass- 
scale, too, suppresses diagrams that violate the OZI rule. 

The chiral behaviour we predicted for the leptonic de- 
cay constant of the pseudoscalar meson excited states is 
exemplified in Fig.|3J viz., f Wl = for m = 0. It is no- 
table that this decay constant is negative as the current- 
quark mass increases away from the chiral limit. Hence 
the residue in the pseudovector vertex of the pole asso- 
ciated with the 2 1 5'o meson is negative on rh > 0. It is 
evident from the table that the residue at the associated 
pole in the pseudoscalar vertex, p Wl (Q , is also negative. 
Thus is Eq. (|16|l satisfied. We anticipate that the pole 



TABLE I: Calculated results for properties of the ground 
and first excited state pseudoscalar mesons: wi(Co) '-~ 
m/(lnCo/A QC D) 7m , Co = 1 GeV; C = 19GeV. Available ex- 
perimental values (in GeV) 0: m no = 0.14 , m ni — 1.3 ±0.1 ; 



m„ = 0.547, m v( 2S) = 1.293±0.005; f„ = 0.092 . Dimensions 
in table are GeV except for /9(C), which is listed in GeV . 



ra(Co) 


m 7T[t 


m 7ri 
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1.08 


0.091 





(0.51) 2 


-(0.49) 2 


0.0055 


0.14 


1.10 


0.093 


-0.002 


(0.52) 2 


-(0.49) 2 


0.125 


0.69 


1.40 


0.130 


-0.023 


(0.64) 2 


-(0.54) 2 



residues alternate in sign; i.e., they are positive for even n 
but negative for odd n. This is required of a vertex that, 
considered as a function of P 2 , is continuous and does not 
vanish between adjacent bound state poles. Since it is the 
square of the residues that appear in the scalar functions 
which characterise the axial- vector and pseudoscalar vac- 
uum polarisations: Ha(P 2 ), Up(P 2 ), this feature is un- 
likely to have readily observable consequences. 

The parameter u> in Eq. I|27|l defines a length-scale 
r a = 1/oj, and magnifying r a increases the range of strong 
attraction in the model. In our calculations we found, as 
anticipated, that the properties of the pion's first radial 
excitation are sensitive to the long-range part of the in- 
teraction. For example, while a 30% increase in r a raises 
m„ by merely 3%, it reduces m 7Tl by 14%. In contrast, 
on the domain of current-quark masses for which our 
present calculations are reliable, the results are consis- 
tent with m no /m 7ri — ► 1 _ in the heavy-quark limit. 

We estimated the charge radius of the pion and its first 



radial excitation using the symmetry preserving impulse 
approximation introduced in Ref. J28J with the Ball-Chiu 
Ansatz [23 for the dressed-quark-photon vertex: r ni = 
1.7 rv . While a more reliable calculation would employ a 
self-consistently calculated drcsscd-quark-photon vertex 
30], our estimate for the ratio should be a fair guide. 

It is natural to extend our quantitative analysis to: 
larger current-quark masses; higher radial excitations; 
and systems in which the constituents have different 
masses. This expansion of the domain on which the man- 
ifestly covariant and symmetry preserving model is ap- 
plied will inform our interpretation of observables, and 
guide the model's improvement and further employment. 
It will not, however, alter the exact results the model has 
been used to illustrate. They serve as a constructive con- 
straint on approaches advocated for the study of (exotic) 
hadron spectroscopy and interactions. 
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